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The curved spacetime surrounding a rotating black hole dramatically alters the structure of 
nearby electromagnetic fields. The Wald field which is an asymptotically uniform magnetic field 
aligned with the angular momentum of the hole provides a convenient starting point to analyze the 
effects of radiative corrections on electrodynamics in curved spacetime. Since the curvature of the 
spacetime is small on the scale of the electron's Compton wavelength, the tools of quantum field 
theory in flat spacetime are reliable and show that a rotating black hole immersed in a magnetic 
field approaching the quantum critical value of Bk — m 2 c 3 /(eh) « 4.4 x 10 13 G ~ 1.3 x 10 -11 cm -1 
is unstable. Specifically, a maximally rotating three-solar-mass black hole immersed in a magnetic 

O ' field of 2.3 x 10 12 G would be a copious producer of electron-positron pairs with a luminosity of 

O ' 3 x 10 52 erg s _1 . 
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I. INTRODUCTION 

The recent discovery that gamma-ray bursts (GRBs) are associated with galaxies at cosmological distances has 
spurred the development of thoeretical models of the central engines of these objects which emit ~ ]_q51-53 er g S 
(assuming isotropic emission) over the span of several to several hundred seconds. Some of the more popular models 
involve some sort of electromagnetic bomb: a quickly rotating, strongly magnetized neutron star (e.g. or a 

rotating black hole threaded by a strong magnetic field (e.g. The first model extends the standard picture of 

radio pulsar spindown (e.g. |6|) to ultrastrong magnetic fields and high spin frequencies; the second model recasts the 
Blandford and Znajek mechanism for the central engine of quasars in the realm of a stellar black hole accreting 
| the debris of a tidally disrupted neutron star. 

An examination of the instability of the magnetized vacuum surrounding a rotating black hole provides an excellent 
starting point to understanding these processes. Van Putten has studied the analogue to Hawking radiation for a 
rotating, magnetized black hole and finds that if the applied field approaches the quantum electrodynamic critical 
^jjy value of 4.4 x 10 13 G a rotating stellar-mass black hole will produce ~ 10 49 erg s^ 1 in pairs. Although this 

technique based on Hawking radiation provides an estimate of the pair production near the hole, the pair production 
for a strongly magnetized, stellar mass black hole depends only extremely weakly on the Hawking temperature of the 
hole or equivalently on the spacetime curvature near the horizon; therefore, accurate results may be obtained if one 
ignores the effects of spacetime curvature on the quantum mechanics of the electromagnetic field surrounding the hole. 
Specifically, Gibbons argues that if the mass, M, of black hole greatly exceeds 10 17 g, the quantum mechanical 
effects of spacetime curvature may safely be ignored for particles more massive than an electron. 

Gibbons examined the problem of how an uncharged rotating black hole embedded in a magnetic field will 
acquire a charge (Q — 2BJ, the Wald charge ]T^] ) through pair-creation near the horizon. The current paper builds 
on Gibbons's picture |ll[| and examines in detail the pair-creation process after the hole has acquired the Wald charge. 
The spacetime curvature does affect the structure of the applied electro magn etic field on scales comparable to that 
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of the black hole. This paper begins with a treatment of this effect in § [I A t hrou gh a discussion of the the Wald 
|l2| field for spacetimes which admit both timelike and axial Killing vectors, jj tIB| specializes this discussion to the 
spacetime surrounding a rotating black hole, the Kerr spacetime |12|-|l4||. § II C calculates the pair production rate 



in a locally inertial frame threaded by both an electric and magnetic field based on the Heisenberg-Eulcr lagrangian 



[|1 5|— |l 7|| - This theoretical basis is utilized to calculate both numerically (§ III A ) and analytically (§ III B ), the pair 
production near a rotating, magnetized black hole. 
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II. THEORETICAL BACKGROUND 



A. The Wald Field 



In a vacuum spacetime, a linear combination of Killing vectors yields a solution for the electromagnetic vector po- 
tential also in vacua. The spacetime surrounding a rotating blackhole yields two Killing vectors ip and rj, corresponding 
to rotations about the angular momentum axis and time translations. Wald found that if the electromagnetic 
field asymptotically becomes a uniform magnetic field, the vector potential is given by 

^ = ^o(Wf^)-^ (1) 

where Q is the charge of the hole. The electromagnetic field is assumed to be a test field, i.e. it does not curve 
spacetime. If the scalar potential of the horizon differs from that at infinity, the hole preferentially accretes charge 
until the potentials are equal. This occurs for Qw = 2BqJ and 

A, L = \b %. (2) 

During the production of the electron-positron jets, the charge of the hole may depart from the Wald value; therefore, 
if Q' = Q — Qwi the complete vector potential is given by 



A n = T^oVv ~ 2jtf V ( 3 ) 



B. The Kerr Geometry 

Since the calculation of the electromagnetic field near the black hole focusses on the Killing vectors, it is propitious 
to use the Boyer-Lindquist coordinates in which the metric takes the form p2|-p^|, 



ds = — I 1 — 
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dt 2 - 



4Mar sin 



-dtdcj) - 
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sin 



^dr 2 
A 



Y,dd 2 



(4) 



where 



E = r 2 
A = r 2 



i 2 cos 2 9 



a z - 2Mr 



(5) 
(6) 



and the Killing vectors are rf = [1, 0, 0, 0] and ^ = [0, 0, 0, 1]. In these coordinates, the field tensor is simply related 
to the derivatives of the metric. If the charge of the hole differs from the Wald value, the tensor consists of two 
components, 



^ 2 ^ 2M » v 



Q' 



(7) 
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303,1 303,2 

-3o3,i 

"303,2 

333,1 333,2 





"333,1 
"333,2 




(8) 



and 
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The invariants / and J also depend simply on the metric coefficients and their derivatives, 
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where g = det(g llv ) and 
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It is important to verify that the various quantities are well behaved on the horizon. A useful expression is the 
value of I at the pole, 



/ = 2B 2 



Q' 



r H 



M 2 



(20) 



J vanishes over the entire horizon and r# = M + \J M 2 + a 2 , the radial coordinate of the horizon. 



C. The Effective Lagrangian of Quantum Electrodynamics 



For a uniform external field the effective Lagrangian of quantum electrodynamics may be written as ^| 
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2 B k 



where -2(ct + ib) 2 = I + I = F^F^ 
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(21) 



T^F^, B k = m 2 c 3 /(eh) w 4.4 x 10 13 G w 1.3 x 1CT 11 cm" 1 and 



is the mass of the electron, a and b are the strengths of the electric and magnetic field measured in a reference frame 
where the two fields are parallel and transform as scalars. 

The pair production probability (w) is simply related to the imaginary part of the Lagrange density of the electro- 
magnetic field p6| , p7[ , w = 2(4ttTi)^ 1 '<sC. The imaginary part of the integrand in Eq. ^l] is even along the real axis 
so the range of integration may be extended over the entire real axis and the contour completed with a semicircle 
encompassing the negative imaginary portion of the complex plane. The integrand has poles along both the real and 
imaginary axes; however, since the integrand is even along the real axis, the residues for the real poles cancel in pairs 
leaving |18|, 
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The more familiar limiting case is where !i«a which yields |T( 
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Taking the opposite limit yields 



h 

mc 



mc* 



^2 dilog [1 



exp(-7nB fe /|a|)] 



(23) 



^log[l 
k 



exp(-7i\B fe /|a| 



(24) 



For a given value of a the pair production rate increases monotonically with b and linearly for b^> a. 
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III. PAIR-PRODUCTION NEAR ROTATING BLACK HOLES 



If the mass, M, of black hole greatly exceeds rnp lanck /m « 10 17 g, the tools of quantum field theory in a flat 
spacetime are adequate to describe the pair production rate near the black hole ; combining the results of the 

previous sections yields a definitive prediction for the pair production and emission energy from the vicinity of the 
black hole before the magnetosphere forms. 

If the magnetic field is parallel (antiparallel) to the angular momentum, positrons (electrons) tend to escape to 
infinity, and the hole quickly acquires a slight negative (positive) charge, so that equal numbers of each charge escape 
to infinity. For a maximally rotating hole, the bulk of the pair creation occurs between latitudes of 30° — 50°. If the 
magnetic field is parallel to the angular momentum of black hole, the positrons escape from the lower half of that 
range. For more slowly spinning holes, the emission region moves closer to the equator. The emission rate on the 
horizon itself vanishes unless the hole has a significant amount of charge (To), i.e. Q' /M 2 < (—4 + 2V3)Bq. 




FIG. 1. The values of the Lorentz invariant "electric" and "magnetic" components of the electromagnetic field surrounding 
a maximally rotating black hole with the Wald charge, Qyj = 2Bq.J. 

Fig. [j] demonstates that the local strength of the electric field inside the static limit is comparable to that of the 
applied magnetic field. In regions where a ~ B^, the pair-production rate is on the order of 10 52 cm -3 s _1 or L ~ 10 2 
(in the rest mass of the particles alone) over the entire spacetime. Clearly, pair production near a rotating black hole 
must become important for a <C Bk 

A. Numerical Results 

Summing the pair production rate over a coarse grid yields an an approximate picture of the pair production near 
rotating black holes. Assuming that the particles reach infinity with the electrostatic injection energy from the region 
where they appear, gives an estimate of the total pair production luminosity. 
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FIG. 2. The pair-production rate and luminosity for a black hole with a = 0.5M and Bo = 0.27Bfc 

Fig. |^ depicts the pair-production rate for a black hole with a = 0.5M with Bq = 0.26£?fc. The charge of the black 
hole deviates from the Wald value with Q' /M 2 = —0.2Bk to make the outflow approximately neutral. The total 
charge of the hole is Q = 0.23Qw- The pair production is highly concentrated near r = 2.1M and a latitude of twenty 
degrees. This is significantly above the horizon which lies at r ss 1.86M and slightly outside the static limit. The 
luminosity is double peaked since the zero potential surface runs through the peak of the pair production. 



TABLE I. Summary of pair-production luminosities from subcritically rotating black holes (numerical results for 
L = 10 50 (M/lkm) 4 erg s _1 ). 
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Since the pair production is highly localized, especially for slowly rotating black holes, the estimates of the total 
luminosity are rather sensitive to the resolution of the grid (higher resolutions may yield higher luminosities), and it 
is difficult to estimate the value of Q/M 2 required to achieve strict charge neutrality in the outflow. 



B. Analytic Treatment 



Since the pairs are produced in a small region of spacetime around the black hole where a <C M, several important 
simplifications are available. First only the first term in Eq. |2^ will be important. Second, the pair-production rate 
near the peak is approximately Gaussian with characteristic widths in the r— and 6— directions. Third, in this small 
region where the pair production peaks, spacetime curvature can be neglected. Fourth, in the vicinity of the peak, 
gradient of the electrostatic potential is constant in magnitude and perpendicular to the zero potential surface, so the 
passage of the zero potential surface through the peak itself guarantees the charge neutrality of the outflow. 

The expansion of the pair production rate about the peak is straightforward and yields, 



"' = Wp 0ak exp AnvAr^Ar 1 ' 



(25) 



where Ar' 1 = 



and A^y is positive semidefinite (its nullspace consists of t — <f> plane) and its other eigenvalues 



are large compared to the characteristic wavenumbers of the blackhole, l/M and 1. The electrostatic injection energy 
relative to infinity vanishes at the peak and it also may be expanded near the peak as e = e- 4l Ar^ 

Consistent with the approximations mentioned earlier, it is also immediate to integrate the pair production rate 
and energy flux over space, taking t — Constant slices, 
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where A = A 1 I 2 and <7\02 denotes the product of the nonzero eigenvalues of A. 

The energy released (or expended) as the particles travel from where they form to infinity is given by 



E\ 



^initial = m 



T)^ - 1 



(28) 



e denotes the component of the energy proportional to the charge of the particle. In Boyer-Lindquist coordinates, 
this component is given by 



Q' M 2 r 
M2 £ 



2A/ 2 rsin 2 6» a 



e vanishes on the conical surfaces 
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(30) 



The component of the energy proportional to the mass of the particle insures that within a thin region bounded by 

(31) 



sin 2 9q I 



r 2 + a 2 cos 2 9q m 



Qwr 



particles of neither charge can escape. Near a stellar-mass black hole, this region is on the order of an electron Compton 
wavelength in thickness; therefore, m may be neglected compared to e, leaving only the electrostatic contribution to 
the energy for consideration. The expansion of e near the peak yields 



r sin 2#o 



7q + a 1 cos-^ fc>o 



(32) 
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The location and width of the peak is found by numerically evaluating the pair production rate for a given value of 
Bq and a variable value of Q' /M 2 . Q' /M 2 is varied until the peak coincides with the zero e surface. This procedure 
works well for a < 0.7M. For larger values of a, the peak becomes too broad and splits in two, so the assumptions 
above are no longer valid. However, since the peak is broad, the direct numerical technique outlined earlier gives 
reliable results. 

Fig. ||depicts the value of Bq to produce a luminosity of 10 50 (M/Mq) 4 erg s^ 1 as a function of a for 0.01 < a < 0.81 
and the mean value of 7 for the primary particles. 
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FIG. 3. The left panel shows the strength of imposed magnetic field to produce a luminosity of 10 50 (M/Mq) 4 erg s _1 for a 
given value of a/M. The solid squares give the numerical results. The solid line in right panel traces the mean value of 7 of 
the primary particles and the dashed line traces the approximate formula Eq. 

The value of Bq is well fit by a power-law such that Bq = B c (10a/M)~ 4 / 5 . The luminosity increases exponentially 
with a and superexponentially with Bq . A convenient fitting formula is 



log 
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(33) 



The location of the pair production and typical escape energies of the pairs also changes with the angular momentum 
of the hole. Fig. |^ depicts the colatitude of the peak and the radius of the peak. A comparison of the right panel of 
Fig. H with the left panel Fig. |] verifies that the typical value of 7 depends rather simply on the colatitude of the peak 
and a through Eq. For small values of a (a < 0.6M), one eigenvector of the matrix A points in the ^-direction, 
the value of 7 is approximately given by 



7 ; 



2(70 ar sin 29q Bq M 
yf-K r 2 + a 2 cos 2 9q B c h/mc 



(34) 



where <jq w 0.055. For larger values of a, the eigenvectors rotate away from the r— and 9— directions, and the 
approximation is poorer. 

The pair production peaks near the equatorial plane for small values of a and moves toward the poles as a increases. 
Meanwhile the peak remains outside the static limit until a w 0.7 M then crosses the static limit and moves toward 
the horizon (c.f. Fig. Fig. ^ illustrates how the analytic treatment becomes unreliable for high values of a. As a 
approaches M, the peak becomes broad, and the assumptions which support the analytic treatment become invalid. 
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FIG. 5. The position and principal axes of the peak. For a = 0.1M, 0.72M and 0.81M. The peak crosses the static limit at 
a « 0.727M. 

For a > 0.735M, the value of I on the horizon near the spin axis becomes negative since Q'/M 2 < (-4 + 2\/3)B ; 
consequently, for larger values of a pairs are produced both in the main peak and in the polar regions near the horizon. 
This polar component does not contribute significantly to the total pair production until a ~ 0.8M Fig. || depicts the 
net pair production rate near a black hole with a = 0.8M. The main peak is quite broad with positrons escaping nearer 
to the equator and electrons at higher latitudes. The horizon component of the pair production consists exclusively 
of electrons. Although locally the rate is quite large, only a small volume is active, so its total contribution is small. 
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FIG. 6. The left panel shows the net pair production rate for a = 0.8M with B = 0.165B C and Q'/M 2 = -0.211B C . The 
right panel shows the pair-production luminosity. 

IV. CONCLUSIONS 

The vacuum surrounding a magnetized, rotating black hole is unstable to pair production if the imposed field 
approaches B^ — m 2 c 3 /(eti) » 4.4 x 10 13 G w 1.3 x 10~ n cm -1 . If the mass of the black hole M is much less then 
10 11 cm or 10 6 M , the field does not contribute significantly to spacetime curvature. For the process as outlined to 
operate, the vicinity of the black hole must initially be free of charge, so the large potential gap and strong electric field 
remain stable for subcritical fields; therefore, one would not expect this simple description to apply astrophysically. 

Adapting the model of Goldreich and Julian |J to the case of a rotating black hole provides an estimate of the 
charge density necessary to short the electric field. Two natural definitions of the angular velocity are that of the 
zero-angular-momentum frame at the horizon and at the peak of the pair production. Both choices result in the 
constraint that the electron density must exceed 2 — 4 x 10 15 cm -3 for a one-solar-mass black hole to short the electric 
field. This is several orders of magnitude larger than the Goldreich-Julian density typical for rotating neutron stars. 
Even if the initial charge density does exceed this value, the vacuum case provides important insights. 

Van Putten's estimates for the pair production luminosity fall short of those calculated here by several orders of 
magnitude |^,^| and depend differently on the mass of the black hole and the strength of the magnetic field imposed. 
Since for weak fields the luminosity depends superexponentially on the imposed magnetic field, agreement can be 
achieved in the gross properties of the models by slightly varying the magnetic field strength. For supercritical fields 
the pair production rate locally increases as B 2 , and the electrostatic injection energy, e, is proportional to BqM; 
consequently, the total luminosity increases as BqM 4 - van Putten argues that the total luminosity increases as 
BqM 2 . Reconciling these differences is difficult as this report and van Putten's work treat the underlying physical 
processes differently. 

Some natural extensions to this work are a treatment of the back reaction of the outflow on the spin of the black 
hole and the source of the external magnetic field. By including a more realistic description of the magnetic field far 
from the hole which would likely include a model for its source, the beaming of the relativistic jet could be determined. 
These developments would constrain the duration of the emission, the nature of its onset and its possible modulation. 
The pairs will naturally produce secondary particles as they travel along the curved magnetic field lines. These 
secondaries, the primaries or previously present material could form a magnetosphere around the black hole, causing 
the electromagnetic field to evolve toward a force-free configuration (e.g. p3| ). 

Rotating black holes coupled to strong magnetic fields naturally produce a highly relativisitic, columnated outflow 
with a total luminosity which can easily exceed 10 50 erg/s. The electrons and positrons are initially separated. 
For a > 0.8M, if the applied field is parallel (antiparallel) to the angular momentum of the black hole, electrons 
(positrons) will escape from polar regions near the horizon from a region near the equatorial plane within the static 
limit. Positrons (electrons) escape only from the equatorial region. For small values of a, both electrons and positrons 
escape from a small region outside the static limit whose latitude depends on the value of a. 
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